False Vacuum in the Supersymmetric Mass Varying Neutrinos Model 
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We present detailed analyses of the vacuum structure of the scalar potential in a supersymmetric 
Mass Varying Neutrinos model. The observed dark energy density is identified with false vacuum 
energy and the dark energy scale of order (10~'^eV)'* is understood by gravitationally suppressed 
supersymmetry breaking scale, -F(TeV)'^/Mpi, in the model. The vacuum expectation values of 
sneutrinos should be tiny in order that the model works. Some decay processes of superparticles 
into acceleron and sterile neutrino are also discussed in the model. 
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I. INTRODUCTION 

In recent years, many cosmological observations sug- 
gest that the universe is flat and its expansion is 
accelerating[l|. The most common explanation for the 
question of the origin of such an acceleration is to as- 
sume the existence of unknown dark energy component. 

A dynamical model, which relates a dark energy scalar 
with neutrinos, so called Mass Varying Neutrinos (Ma- 
VaNs) scenario, was proposed as a candidate for the dark 
energy 0, In the context of this scenario there are a 
lot of works 0, H, 0|. Among these works, we consider a 
MaVaNs model in the supersymmetric theory. 

The MaVaNs scenario assumes an unknown scalar 
field so called "acceleron" and a sterile neutrino. They 
are gauge singlets under the gauge group of the stan- 
dard model. In the extension of the scenario to the 
supersymmetry, these fields are embedded into chiral 
multiplets[^, [y]- Typical time evolutions of the neu- 
trino mass and the equation of state parameter of the 
dark energy are presented in Ref. @ , where the acceleron 
and a sterile neutrino are embedded into the single chi- 
ral multiplet. In this model, the observed dark energy 
density is identified with the false vacuum energy and 
the dark energy scale of order (10~^eV)^ is understood 
by gravitationally suppressed supersymmetry breaking 
scale, F(TeV)^/Afpi. These works are based on assump- 
tions that the imaginary part of the acceleron and vac- 
uum expectation values of sneutrinos vanish. 

In this paper, we present a modified supersymmetric 
MaVaNs model of Ref. [6] and detailed analyses of the 
vacuum structure of the scalar potential taking account 
of the finite imaginary part of the acceleron. Since the 
vanishing imaginary part of the acceleron should be dy- 
namically realized, discussions based on the assumption 
in Ref. are inadequate to analyse the vacuum structure 



of the potential. Actually, we find the global minimum 
is different from one in the case of the vanishing imagi- 
nary part of the acceleron. The dark energy is identified 
with the false vacuum energy and the dark energy scale is 
given by pj/^ i^(TeV)VMpi in Ref.0. We show that 
the identification and the observed dark energy scale are 
also realized in the modified model where the acceleron 
is appropriately stabilized. Furthermore, we present cos- 
mological discussions for the magnitude of the Yukawa 
coupling among the acceleron and sterile neutrinos which 
is assumed to be of order one in Ref. [6,]. In this frame- 
work, the model with non- vanishing vacuum expectation 
values of sneutrinos are also discussed. We find that the 
magnitude of VEVs of sneutrinos are strictly constrained 
in the supersymmetric MaVaNs model. Then, two body 
decay processes of the superparticles are also studied. 

In section 2, we present a supersymmetric MaVaNs 
model and detailed analyses of the vacuum structure of 
the scalar potential with finite imaginary part of the ac- 
celeron and vanishing VEVs of sneutrinos. In section 3, 
we discuss the case of non- vanishing VEVs of sneutrinos 
in the model and show constraints on the magnitude of 
VEVs. In section 4, two body decays of the superparti- 
cles into acceleron and sterile neutrino are investigated. 
The section 5 is devoted to the summary. 



II. VACUUM STRUCTURE OF THE SCALAR 
POTENTIAL 

A. Scalar potential for the acceleron 

Let us discuss vacuum structure in a supersymmetric 
mass varying neutrinos model. We consider the following 
superpotential W in the MaVaNs model. 
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where L and R are chiral superfields corresponding to 
the left-handed lepton doublet and the right-handed neu- 
trino, respectively, [m The A is a gauge singlet chiral 
superfield. The scalar and spinor components of A are 
represented by (f> and ip^ respectively. The scalar com- 
ponent corresponds to the acceleron which causes the 
present cosmic acceleration. The spinor component is 
a sterile neutrino. Coupling constants and mass param- 
eters are represented by = 1,2) and M^, M^i, Mr, 
niD , respectively. The first and second terms of the right- 
hand side in Eq. ([1]) are derived from Yukawa interac- 
tions such as yiHLA and y2HLR with yi{H) — mu and 
y2{H) = Md, respectively, where H is the usual Higgs 
doublet. In general, other interactions among gauge sin- 
glet fields such as X'^R^ and X^AR^ are allowed but con- 
tributions from these interactions to the acceleron poten- 
tial and neutrino masses are negligibly small. The terms 
appearing in Eq. ([1]) are important to realize the Ma- 
VaNs scenario. In this model, there is no assignment to 
conserve the R-parity due to the emergence of trilinear 
coupling. In the limit of ?/i — > 0, a dark energy sector 
composed of A is decoupled from the MSSM (visible) sec- 
tor and then R-parity conservation is restored as in the 
MSSM. The LSP is absolutely stable in the case of con- 
served R-parity, however some decays of sparticles into 
the acceleron and a sterile neutrino are generally allowed 
in the model given by Eq. We will return to this 

point later. 

Taking supersymmetry breaking effects into account, 
the scalar potential for the acceleron is given by 
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where m and Ki{i = 1^3) are supersymmetry breaking 
parameters and Vq is a constant determined by imposing 
the condition that the cosmological constant is vanishing 
at the true minimum of the acceleron potential. 

In the Mass Varying Neutrinos scenario, the dark en- 
ergy is supposed to be composed of the neutrinos and the 
scalar potential for the acceleron: 



Pde = Pi/ + ^"(0), 



(3) 



where pde and pi, correspond to the energy densities of 
the dark energy and the neutrinos, respectively. For con- 
venience later, we rewrite the scalar potential ^ in terms 
of two real scalar fields 0^ (i = 1, 2) instead of 0, 



xl + xi 



-^Ma + 2«i j 0102 + ^{cl)l+ 4^1) 
+ (2Mavl + '2MdI'r + A2^>L^^i?j"^D0l 

+{i?L{XimD + X2MD) + X2MRVR}{(j)\ - 4>l) 
-2{k2Vl + K3i>R)H(f>i + Vo, (4) 

where we define as cj) = (pi + i(j>2 and to^ = 2(M^ -I- m|, — 

TO^). 

First we discuss the vacuum structure of this scalar 
potential in the case of the vanishing vacuum expectation 
values of sncutrinos. 
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If {XiMa/2 — 2ki/3) < 0, the scalar potential has a local 
minimum at the origin in field space. Moreover, when 
parameters satisfy a relation, 
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the potential has a second local minimum at 0i > and 
02 =0, V2, which is lower than the local minimum at the 
origin, Vi. Then, V2 is given as follows: 
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In a region of 
_ {XiMa + 4:Ki) + M 

xi + xi 



< 01 < x9 . X9 ,(8) 
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has non-vanishing VEVs, 



, , ^ (A2 + Ai)02 + (AiM^ + 4«)0i+m2 
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where M = ^{XiMa + 4k)2 - 4(A? + Xfiml. The po- 
tential energy of the local minimum in this region, V3, is 
given by 



v. 
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Since V3 is smaller than V2 in the case of XiMa <C ki, V3 
is the value of global minimum of this scalar potential. 
Therefore, we find 
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FIG. 1: The acceleron potential given in Eq. ([5]). There are 
four local minima, the first one is the origin in field space, Vi, 
the second one is situated at > and (j)2 = 0, V2, and the 
remaining two local minima are degenerated and global ones, 
V3. 



In summary, there are four local minima in the scalar 
potential ([5]). They are shown in Fig[T] The first one is 
the origin in field space and the second one is situated at 
> and (/)2 = 0. The remaining two local minima are 
degenerated and global ones of the scalar potential. 



- V{cl>2 - ±{cf>2)g), 
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where 
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We find that a constant Vo in Eq. ([2]), which cancels the 
cosmological constant, is tuned such as Vb = jVaj. 



B. Numerical analysis 

Now, let us consider a possible Mass Varying Neutrinos 
model with this scalar potential. Since the dark energy 
in the MaVaNs scenario is the sum of the scalar poten- 
tial and the energy density of neutrinos, the stationary 
condition is given by 
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We find the local minimum at the origin is slightly devi- 
ated due to the contribution of the neutrino energy den- 
sity. And a numerical calculation shows that the model 
works near the origin in field space: (^i)" — eV and 



k'^'i)^ — eV, where superscript zero means the value at 
the present epoch. Since ^ field is trapped in the local 
minimum near the origin, one observes the potential en- 
ergy IV3I as the dark energy at the present epoch. The 
(\) field trapped near the origin is plausible because it is 
likely that the universe starts at the origin due to the 
thermal effects. Typical values of parameters are given 
as 



Ki ~ 2.17 X lO^^eV, Ma 
\MI/Mr 



1.003 X lO^^eV, 



4.998 X lO^^eV, 
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1.008 X lO^^eV, 
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by numerical calculation when we take 
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as inputs masses and couplings. Now, we assume that 
supersymmetry is broken at TeV scale in a hidden sector 
and its effects are transmitted to a dark energy sector, 
which includes a chiral supcrficld only through grav- 
itational interaction. Then, the scale of soft terms of 
order (TeV)VMpi - ©(lO'^-lO-^ eV) is expected. Soft 
parameters K2 and K3 are three-scalars couplings such 
as K2Vi^li(^ and K^vuHcj). If the dark energy sector is 
communicated with the visible sector, which includes the 
right-handed sneutrino, only through the gravitational 
interaction except for a peculiar Yukawa interaction of 
the MaVaNs scenario such as yiHuLtp^ K2,3 are roughly 
estimated of order m1^^^/M-p\ by dimensional analysis, 
where Wsoft means the scale of soft terms in the visible 
sector. The magnitude of soft parameters in the visible 
sector depends on a mediation mechanism between the 
visible and a hidden sector. If we take rn-soft ©(TeV), 
the magnitude of K2.3, becomes 0(10"^ eV). A larger 
scale of K2,3, which corresponds to msoft, prefers tiny or 
vanishing VEVs of sneutrinos in the MaVaNs model. We 
will return to this point later. 

The parameter Ma affects on the neutrino masses. The 
effective neutrino mass matrix is given by 



M 



c + ei mD + e2 

m_D + £2 Ma + XiiM - 



£3 



(17) 



in the basis of {vli V'); where i^l and ip are the left-handed 
and a sterile neutrino, respectively. The c in the (1, 1) el- 
ement of this matrix corresponds to the usual term given 
by the seesaw mechanism in the absence of the acceleron, 
c = ~Ml)/Mii, where Ai((/)i) <C Mj^ -C Mr is assumed. 
The parameters £^(1 = 1 ~ 3) are functions of (jji, which 
are suppressed by large scales such as the right-handed 
neutrino mass scale Mr or the higgsino mass . If the 
VEV of (j>i and the mixing between the active and a ster- 
ile neutrino are small, the mass of a sterile neutrino is 
almost determined by Ma- We assume that the mag- 
nitude of a sterile neutrino mass is the same order as 
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the left-handed neutrino one, = 10 ^ eV, and notice 
TO^ ~ Ma ~ 1.003 X 10"2 eV as shown in Eqs. (US]) and 
(USD. 

The parameter mo determines the mixing between the 
active and sterile neutrinos. When we take as m^i = 
10~'^eV, the mixing angle, sin^ ~ 4 x lO^'', is expected. 
The active and sterile neutrinos mix maximally in the 
case of 10~^eV ^ m^, which is ruled out by the current 
astrophysical, cosmological and laboratory bounds, [l^ 
Since rajj is defined as yi {H) in our model, rriD < 10~^eV 
and (H) ~ lO'^ GeV require yi < 10^^*. In general, Ma- 
VaNs models need such a tiny Yukawa coupling among 
the left-handed neutrino, the Higgs boson and the sterile 
neutrino. In the limit of yi — > 0, a dark energy sec- 
tor is decoupled from the visible sector and the mixing 
between the left-handed and sterile neutrinos vanishes. 
Such a small Yukawa coupling may be explained in a 
brane model however the discussion is beyond the scope 
of this paper. We use the value m^i = 10~"^eV, which 
corresponds to yi ~ 10"^^, in numerical analyses. 

In the MaVaNs scenario, the mass of the scalar field 
should be less than O(10"^eV) in order that the scalar 
field does not vary significantly on distance of intcr- 



-1/3 



Therefore, 



m0 



= lO-^eV is 



neutrino spacing, rii, 
taken in our calculation. 

Next, we comment on couplings A^. The size of Ai 
is constrained by the successful Big Bang Nucleosyn- 
thesis (BBN) theory and supernova data. In the early 
universe, the acceleron and the sterile neutrino have to 
decouple in order that the scenario of BBN is not radi- 
cally changed. If the acceleron production rates through 
acceleron-strahlung and pair annihilation of neutrinos are 
less than the expansion rate at the era of BBN, these 
particles are not in thermal equilibrium at the epoch 
These considerations give us 



drrii, 



< 1, 



(18) 



(^) 



Eqs. p8)) and (|19p show constraints that (/)-production 
rates through (/)-strahlung and neutrino pair annihilation 
are smaller than the expansion rate at the BBN epoch, 
respectively. 

Similar constraints come from consideration in super- 
novae. We require that 0-emission does not cool a pro- 
toneutron star too quickly to change the observed neu- 
trino spectra emitted in the first ten seconds. This re- 
quirement is realized in the following conditions 



dull, 



'dm. 



(20) 



< 10 



-23 



'TO. 



< 10-^^eV"^ (21) 



and ((2T|) . are more severe than the ones from the BBN 
scenario, Eqs. pB]) and P^ . 

In our model, the acceleron dependences of masses of 
the left-handed and sterile neutrinos are given by diago- 
nalizing the effective neutrino mass matrix (|17p . 



c + Ma + \i{4>i) 



2 

c + AfA + Ai(0i) 
2 

v/[c-(Ma + Ai(<^i))]2+4to2, 



, (22) 



(23) 



where we ignore e^. Fig [2] shows that the ac- 
celeron dependences of neutrino masses, (dm^/dcfii)^ and 
{d'^m^/d(t>\)'^ , and the constraint on Ai from supernovae. 
Since the neutrino masses depend on only 0i, {d/d(j)) is 
replaced with {d/d(j)i). The left and right figures in Fig[5] 
corresponds to the constraint in Eq. (|20p and the second 
one in Eq. {d^m^/d(p^)'^ < 10~37 eV^^. We find 

that both constraints are satisfied if Ai < 4 x 10^^. 

Next, we consider the coupling A2. In numerical cal- 
culations, we take the same number density of ster- 
ile neutrinos as the one of the left-handed neutrinos, 
rii, oc , where T is the neutrino temperature. This 
means that sterile neutrinos should be thermalizcd in the 
early universe. The interaction rate of a sterile neutrino 
through the interaction (X2/2) ARR is typically given 
by r ~ A|T. The thermalization of sterile neutrinos, 
r > T^/Mpi), requires A^Mpi > T. For instance, 

if we assume that sterile neutrinos thermalizcd after the 
inflation, X^Mpi > Tjih should be realized, where T^h is 
the reheating temperature. The thermalization of sterile 
neutrinos in this MaVaNs model with the coupling A2 of 
order one is easily realized by many inflation models. 



C. False vacuum 

The presented MaVaNs model works at the false vac- 
uum. We estimate the decay probability of the false vac- 
uum. The decay probability is given by , where B is 
the action for the classical solution in Euclidean space. 
We use the result of 0], 



B = 2n^ 



(A0)' 



[(Ay_)i/3 _ (Ay+)i/3]3' 



(24) 



where Ac/) = 0+ , 0+ and 0_ are the vacuum expec- 
tation values at the false and the true vacuum, respec- 
tively, and 



AV±=VT~Viq}±), 



(25) 



We find that constraints from supernova data, Eqs. (I^Dl) 



where Vr is the potential energy of the top of the barrier 
between the false and the true vacuum. After calculating 
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FIG. 2: The acceleron dependences of neutrino masses and 
the constraint on Ai from supernova data. The left figure 
shows the relation between [dm^ / d(f)Y and Ai. The horizon- 
tal line shows the constraint of Eq (|20p . The dashed and 
thin curves correspond to (9m^^/9(^i)^ and {dm^ / dcjii)^ , 
respectively. The lower region than the horizontal line, 
{dm^/d(l)if < 10"^^ is allowed by the supernova data. In 
the right figure, the horizontal line shows the constraint of Eq. 
(|21|) . The thin curve corresponds to (d'^m^j^/dtfii)^ , which is 
exactly equal to {d^m^ /dcj>i)^ . The lower region than the 
horizontal line, {d^Tn^/d(f>^)^ < 10"'^^ eV~^, is allowed by 
the supernova data. These figures mean that if Ai is smaller 
than about 4 x 10~^, all constraints are satisfied. 



the dark energy density numerically in our model, we find 
6+ ~ 3.61 X IQ-^eV, (h- ~ 2.83 x IQ-^eV, 



AV+ ~ 2.3 x 10"2"eV^, AV- ~ 2.95 x 10~"eV^, 

(26) 

and thus B ~ 42.9 ^ 1. Therefore, the decay probability 
is sufficiently small. This probability is of the decay from 
the false vacuum near the origin in the field space, Vi, 
into one of two degenerated true vacua, V3. There is 
another possibility of the false vacuum decay, which is 
the one into another false vacuum, V2- In this case, we 
get B ~ 41.1, where we take 0_ ~ 2.80 x 10"^ eV. In 
this numerical calculation, we take the square potential 



approximation in Ref. \f\. In the limit of AV+ 0, the 
result of this approximation Eq. ([M)) describes the sheer 
drop which is the extreme case of "tunnelling without 
barriers" discussed in Ref. 



B = 27r 



2 iMr 

AV- ' 



(27) 



Since AV4. <§; AV- in our model, this approximation 
is appropriate rather than "triangle" ff\ or "thin-wall" 
approximation [9]]. 



III. VACUUM EXPECTATION VALUES OF 
SNEUTRINOS 

In the previous section, we consider the vanishing vac- 
uum expectation values of sncutrinos. Let us consider 
effects of non- vanishing VEVs of sneutrinos. VEVs of 
sneutrinos interacting with the acceleron are strictly con- 
strained in the Mass Varying Neutrinos model. The most 
severe one comes from the stationary condition Eq. ^ 
We rewrite the condition Eq. as 



dV 
d4>i 



.^3 ^ 

I = llZ^,Tp 



dmi dFj^i) 



(28) 



where = T^J^Fi^i), 6 = mi/T and 



1 r 



ev + 1 



(29) 



Since the right hand side in Eq. (^5]) is approximately 
proportional to the number density of neutrinos and both 
active and sterile neutrino masses derived from the mass 
matrix Eq. p7p do not depend on (j)2, the stationary 
condition at the present epoch is approximated by 



dV 



dmi 



(30) 



01 = (01>O. 



A viable MaVaNs model should satisfy this constraint. 
Of course, the model with the vanishing VEVs of sneu- 
trinos satisfies it. When we take Ai — 10~®, we have 



0(10"'^), 



0(10" 



(31) 
(32) 



Therefore, the magnitude of the right hand side in Eq. 

©(IQ-i^eV^). 



30]) is O(lO-^^eV^), where we take 



Since the term 



61 is dominant in the left hand side of 



Eq. (1201), we have (01)° ~ O(10-"eV), where we take 
= 10^'' eV. 

When the vacuum expectation values of sneutri- 
nos are non-vanishing, the following terms are strictly 
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constrained unless miraculous cancellation among each 
terms is realized, 



+XiX2VR(l)l + 2\2Mavr4'i 
+2{£>L(Aimi) + \2Md) + \2MRVR}(j)i 
~2{k2Vl + >i?.VR)H. (33) 

Especially, restrictions on the VEVs of sneutrinos are 
obtained from k.2Hvl and n^iHuR. In order that these 
terms do not spoil the model, {vl.r) < 10^'^^ eV must be 
realized when we take ^2,3 = 10~'^eV and (H) = lO^GeV. 
The larger values of ^2,3 gives the smaller VEVs of sneu- 
trinos. The vanishing or tiny VEVs of sneutrinos are 
favored in this supersymmetric Mass Varying Neutrinos 
model. 



IV. SPARTICLE DECAY 

The Mass Varying Neutrinos scenario assumes some 
light particles such as the acceleron and the sterile neu- 
trino. In the supersymmetric realization as we discussed, 
these particles are embedded into single chiral multiplet. 
Since Eq. ^ leads to R-parity violation in order that 
a given scalar potential and the neutrinos mass are con- 
sistent with the current cosmological observations, some 
heavy sparticle such as the sneutrino and the higgsino 
decay into light particles due to such R-parity violating 
interactions. 

We consider two body decay processes of the sneutrino 
and the higgsino in this section. A peculiar decay process 
of the sneutrino in the MaVaNs scenario is i>L — * 
whose decay rate is given by 



D 



(34) 



When we take mp^^ — ITeV, Ai = 10~ and m^) ~ 
lO^^eV, we have l/r£>^^00 ~ 3 x lO^^sec. A larger left- 
handed sneutrino mass or smaller Ai and mo suppress 
this process. The dependence of on the left- 

handed sneutrino mass and Ai is shown in FiglS) 

The higgsino decay is another peculiar one in this sce- 
nario. There are two important processes, h — > VLcf) and 
The decay rates of these processes are given by 



327r ' 



r 



327rm? 

h 



(35) 



ITeV, l/r,;_^ 



and m,~ - ...v, 
If we assume that the sneutrino mass is less 

is 



When we take yi — 10 
700 sec 

than the higgsino mass but the same size, l/E^^^^^^ 
the same order as l/E^^^^^^^. Larger uif^ or yi enhances 
these processes as shown in Figl3l The above processes 
as we discussed may affect on many LSP cold dark mat- 
ter models. Such effects of the Mass Varying Neutrinos 
scenario will be discussed elsewhere. 
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FIG. 3: The left figure shows the dependence of l/Fpj^^^^ 
on the left-handed sneutrino mass and Ai, where mo = 10""^ 
eV is taken. The region lower than the bold horizontal line, 
which corresponds Ai ~ 4 x 10~^, is allowed. The right figure 
shows ~ l/Tf^^^^^ ~ yfmf^/32TT. 



V. SUMMARY 

We presented a supersymmetric Mass Varying Neutri- 
nos model and detailed analyses of the vacuum structure 
of the scalar potential taking account of the finite imag- 
inary part of the acceleron. The potential has four local 
minima and the dark energy is identified with the false 
vacuum energy. A metastable vacuum are realized near 
the origin in field space, ~ and <p2 — 0, due to the 
supersymmetry breaking effect. This metastable vacuum 
is enough stable in the age of the universe. In the model, 
the observed dark energy scale of order (10"^eV)^ is un- 
derstood by gravitationally suppressed supersymmetry 
breaking scale, i^(TeV)V^Pi- 

Following these considerations, the case with non- 
vanishing vacuum expectation values of sneutrinos have 
been discussed. We have found that the VEVs of sneu- 
trinos are strictly constrained such as {vl,r) < I0~^^ 
eV in order that terms with respect to sneutrinos in the 
acceleron potential do not spoil the model. This means 
that the vanishing or tiny VEVs of sneutrinos are fa- 
vored in the proposed supersymmetric MaVaNs model. 
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Finally, two body decay processes vl 4>4'i h ^ v^cf) 
and h ^ i^i^tp are discussed. The obtained decay rates are 
^ 3 X 1022 sec. and l/F^^,^^ c l/r^_,,^^ ^ 
700 sec. Effects on many LSP cold dark matter model 
from these processes is expected to provide interesting 
physics. 
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